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Making chiral fermion actions (almost) gauge invariant using Laplacian 
gauge fixing 
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a University of California at San Diego, Department of Physics, La Jolla, CA 92093, USA 

Straight foreward lattice descriptions of chiral fermions lead to actions that break gauge invariance. I describe 
a method to make such actions gauge invariant (up to global gauge transformations) with the aid of gauge fixing. 
To make this prescription unambiguous, Laplacian gauge fixing is used, which is free from Gribov ambiguities. 



1. LATTICE CHIRAL FERMIONS 

There are many proposals to describe chiral 
fermions on the lattice (see refs. [0J|] for re- 
cent reviews), which all appear to work when 
the fermions couple to smooth external gauge 
fields. When the gauge fields are dynamical, how- 
ever, these approaches fail: e.g. because oppo- 
site handedness mirror fermions are generated dy- 
namically H, or because the chiral fermions are 
shielded from the gauge fields and no longer cou- 
ple to it |^|Jj|. The lesson learned from these 
attempts is that the 'longitudinal mode' of the 
gauge field, which correspond to gauge transfor- 
mations, is responsible for the breakdown of these 
models and must be controlled. 

To constrain the dangerous longitudinal mode, 
the Rome group has proposed to include gauge 
fixing in the definition of the path integral Q. 
Now the dynamical, but gauge fixed gauge fields 
should be sufficiently smooth on the lattice scale 
such that the chiral fermions remain chiral, as 
was the case with classical external gauge fields. 
To recover the gauge invariant target model it 
should then be sufficient to include counter terms 
in the action and tunc their coefficients such that 
BRST invariance is recovered. This approach 
looks promising, if it can be implemented non- 
perturbatively. This is difficult, however, since 
gauge fixing requires a Faddeev-Popov determi- 
nant which is not positive definite and also Gri- 
bov ambiguities may pose problems. 

Here I discuss a practical method to implement 
the Rome idea, without having to put gauge fix- 
ing and ghost terms in the action. As an example 



I consider Wilson fermions with left coupling to 
a gauge field U 9 , 

S = ^ x mU 9 UPL+^yPR)^y (1) 

xy 

X 

wiihp and the gauge covariant and free lattice 
Dirac operators respectively, Pr : l — (1 ±7s)/2. 
The Wilson term breaks the left-gauge invariance, 
ipx — ► (Pr + 9xPh)i>x, g x G G. However, the 
breaking term is oc a and should be irrelevant for 
a — > 0, at least for classical fields. 

The link field U 9 is a gauge transform of U. 
When integrating over all U in the path integral 
(without gauge fixing) it is equivalent to inte- 
grate with measure DU, DU 9 or DUDg. Then 
the field g represents the longitudinal gauge mode 
and could be interpreted as a spurious Higgs field, 
which should not affect the physics. However, 
this non-gauge fixing approach was found not to 
work • In the Rome approach one adds gauge 
fixing and ghost terms (and counter terms, of 
course) , such that U 9 is constrained on a smooth 
gauge section. 

In the present approach U 9 also satisfies 
a smooth-gauge condition. However, this is 
achieved not by putting a gauge fixing term in 
the action, but by actually computing the gauge 
transformation g(U) which fixes U . In the path 
integral one still integrates over all U, but the 
field to which the fermions couple is not U but 
U g ( u K Assuming that g(U) is uniquely de- 
fined for all U, it follows that under a gauge 
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transformation U — ► U h , g(U) transforms as 
g(U h ) x — h g(U) x h x , with h a global trans- 
formation and one can verify that the action (^) 
is invariant up to a global transformation with 
ho, S(ip,lp,U h ) = SO,?,^' 10 )- Notice that the 
fermion fields are not transforming. With the uni- 
tary (gauge) transformation ip — > (Pr + gP^)^, 
ip — > ip(g'Pn + Pl) the g can be removed from 
Ip and the fermion fields transform now as usual. 
For abelian gauge groups U h ° = U and then the 
action is actually gauge invariant (see ref. j| for 
the details). 

Several remarks are now in order. 

1) It should be stressed once more, that the ac- 
tion (||) with g = g(U), is a function over the full 
gauge field configurations space. The integration 
is with the usual Haar measure DU, not with the 
Faddeev-Popov measure that includes gauge fix- 
ing and ghost terms. Note that counter terms 
should still be included in the action. 

2) The field g(U) depends on all link variables 
Ufj,x-i and therefore the action is nonlocal. I as- 
sume here that at least for anomaly free fermion 
representations, the action effectively becomes lo- 
cal in the scaling region. 

3) The gauge transformed field must be suffi- 
ciently smooth, such that one can hope that the 
good behavior found with classical U fields also 
persists for dynamical fields U 9 ^ u '. 

4) It must be possible to compute the gauge trans- 
formation g(U) without Gribov ambiguities, in a 
reasonable amount of computer time. 

5) For practical purposes it is essential that the 
nonlocal action can be simulated efficiently. 

In ref. || points 4) and 5) are addressed and 
it is shown there that a hybrid Monte Carlo al- 
gorithm can be used to simulate the action. The 
details depend on the choice of gauge condition. 
In order to comply with point 4), the Laplacian 
gauge can be used, which was introduced in ref. 
0. Here g(U) is computed from the eigenfunc- 
tion of the gauge covariant Laplacian with the 
smallest eigenvalue (for G =U(1) or SU(2)). For 
the HMC algorithm it is then required to com- 
pute this lowest eigenfunction numerically and 
to invert a Laplacian-like matrix, at each HMC 
time step. This is time consuming, but when 
the fermions are dynamical anyway, the addi- 



tional work is relatively small. The Laplacian 
gauge is unambiguous, provided that this lowest 
eigenvalue is nondegenerate. A second singular- 
ity arises when the eigenfunction would vanish at 
some lattice site. Even though these 'Gribov hori- 
zons' are of measure zero in the path integral, one 
should check point six: 

6) Gribov horizons should not pose a problem in 
practice. 

2. SOME TESTS 

Points 3) and 6) raised above can be addressed 
to some extent, using quenched dynamical gauge 
fields. This is done in ref. pp] , which contains 
the details. A measure for the smoothness of the 
gauge fixed field is the value of the average link 
(U). The Landau gauge maximizes this quan- 
tity. It turns out that the results of the Landau 
and Laplacian gauges are very similar. Fig. 1 
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Figure 1. Average link after standard Lan- 
dau (solid), Laplacian (dashed) and Laplacian fol- 
lowed by Landau gauge fixing (dotted), for U(l) 
gauge fields on a 20 2 lattice. 

shows the average link in the U(l) model in two 
dimensions. It is seen that Gribov copies (local 
maxima) of the Landau gauge make (U) actu- 
ally slightly smaller than in the Laplacian gauge 
for j3 > 3. Similar measurements with dynam- 
ical SU(2) gauge fields in four dimensions show 
even less difference between (U) in the Landau 
and Laplacian gauges. 
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Even if the Laplacian gauge leads to equally 
smooth gauge fixed fields as the Landau gauge, 
one can still fear that this is not smooth enough 
for the fermions. In particular when there is a 
gauge singularity, as in the core of an instanton, 
the U 9 ^ field locally has to deviate far for one 
and this might give problems for the fermions. 
To test for this, we have computed the eigenvalue 
spectrum of the fermion matrix defined implicitly 
in eq. @ , using test configurations with nonzero 
topological charge. With these backgrounds the 
fermion matrix should have a zero mode, cf. ^J§. 
It is found that the left-coupled Wilson fermion 
matrix has an (almost) real eigenvalue, with a 
small nonzero real part which is found to scale 
to zero oc a. This shows that in the continuum 
limit, the fermion matrix indeed has its desired 
zero mode. 
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Figure 2. Flow of the smallest two eigenvalues in 
the 2d U(l) model at (i — 4 on an 8 x 20 lattice. 

Also point 6) is addressed in ref. jlQj . Using 
the U(l) model in two dimensions, I have com- 
puted the lowest two eigenvalues of the Laplacian 
after each time step of a HMC simulation of the 
quenched model. For values of the gauge cou- 
pling in the scaling region, /3^3, the separation 
between the lowest two eigenvalues only becomes 
comparable with the fluctuations of the individ- 
ual eigenvalues when the lattice is (unnecessarily) 
large in at least one direction. It is found that the 
eigenvalues occasionally approach each other, but 
never come so close that the computer code runs 



into trouble computing them. As an example we 
show part of the level flow found at f3 — 4 on an 
8 x 20 lattice. 

Also in SU(2), where I monitored the level flow 
during Metropolis updating, I only found avoided 
crossing, with a smallest separation between the 
eigenvalues which is many orders of magnitude 
larger than the precision with which these eigen- 
values can be computed. Also the other kind of 
singularity, the vanishing of the eigenfunction at 
some site, never came close in these test runs. 
These results indicate that the Laplacian gauge 
can be computed very efficiently and without nu- 
merical problems. Hence it seems worthwhile to 
further investigate this approach to chiral models 
on the lattice 
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